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The goal of this research is to accurately predict the characteristics of the laminar 
separation bubble and its effects on airfoil performance. To this end, a model of the 
bubble is under development and will be incorporated in the analysis section of the Eppler 
and Somers program. 1 As a first step in this direction, an existing bubble model has been 
inserted into the program. This has the dual purpose of gaining insight into the correlation 
between existing theory and bubble prediction trends, as well as becoming familiar with 
the analytical and numerical details of Eppler’s boundary-layer analysis method. 

The selection of the first bubble model to be used in the Eppler and Somers program 
was based on several considerations. In the first place, it was decided to address the 
problem of the short bubble before attempting the prediction of the long bubble. The 
short bubble entails, by definition, a much weaker interaction between the viscous and 
inviscid regions of the flow field than the long bubble. This, and the greater amount of 
experimental data available for the short bubble, renders the development of a short-bubble 
model more tractable both analytically and numerically. In the second place, an integral 
boundary-layer method is believed more desirable than a finite difference approach. While 
these two methods achieve similar prediction accuracy, finite-difference methods tend to 
involve significantly longer computer run times than the integral methods. Since the Eppler 
and Somers program is primarily an airfoil design program, computational efficiency of 
the bubble model is of primary importance. Finally, as the boundary-layer analysis in the 
Eppler and Somers program employs the momentum and kinetic energy integral equations, 
a short-bubble model compatible with these equations is most preferable. 

* NASA Technical Monitor for this Grant is Mr. Dan M. Somers, NASA Langley 
Research Center, Hampton, VA 23665 
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The short-bubble model of Horton 2 meets with all the above considerations and was 
therefore chosen for incorporation into the program. In addition, some modifications based 
on experimental data more recent than that available to Horton have been included. Specif- 
ically, the length of the laminar length is based on a curve fit, calculated by O’Meara and 
Mueller, 3 of Reynolds number based on the laminar length as a function of Reynolds num- 
ber based on the momentum thickness at separation. Values suggested by Roberts 4 for the 
mean dissipation coefficient in the turbulent separated shear layer and for the reattachment 
parameter have been used. The bubble is matched to Eppler’s laminar separation crite- 
rion and turbulent boundary-layer method. In spite of these improved correlations and 
empirical constants, this model remains quite crude. As will be shown below, however, 
this crudeness has helped to isolate those areas in the bubble region and in the integral 
boundary-layer method that will need further study before a faithful reproduction of the 
pertinent physical processes can be achieved. 

Analytical and Numerical Details 

The modified Horton bubble model which has been incorporated into the Eppler and 
Somers program is triggered when Eppler’s laminar separation criterion is met; that is, 
when 

(jffsi)s = 1.515095 (1) 

As discussed by Horton, in the laminar portion of the bubble the streamwise pressure 
gradient is assumed to equal zero. Thus, 

U T = U S (2) 

Further, the assumption of negligible skin friction leads, from the momentum integral 
equation, to 

(S 2 )t = (fo )s (3) 

O’Meara and Mueller have correlated the length of the laminar shear layer to the momen- 
tum thickness at separation, based on seven different sets of data including their own, and 
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found 


( 4 ) 


Re e i = 155/Ec(5 2 ) 5 


or, 

h. = 155 <M£ 

c c 

This correlation matches well more recent data taken by Brendel and Mueller, 5 and was 
therefore used in place of Horton’s original transition criterion, 

Re ti = 40000 (5) 


In order to proceed from the transition point to the turbulent reattachment point, 
a number of additional assumptions are made. As discussed in Ref. 2, it appears from 
experimental data that 


dHz 2 

ds 


= 0 


( 6 ) 


Based on this result, and on the vanishing of skin friction at a point of reattachment, 
Horton was able to reduce Truckenbrodt’s shape parameter equation to 


'< 5 2 dU' 


Cd 

Uds 

R 

[#32(#12-1)J 


= A R 


( 7 ) 


The assumption of a universal reattachment velocity profile and a constant Cp in the 
turbulent part of the bubble leads to a constant value for Ar. Roberts suggests an exper- 
imentally determined mean value of the dissipation coefficient, 


C D = 0.035 


( 8 ) 


Using the generally accepted values of 


{H Z2 )r = 1.51 


( 9 ) 


and 


(# 12 )/? = 3.5 


( 10 ) 


3 



Eq. (7) gives 


A R = -0.0093 


( 11 ) 


As discussed in Ref. 6, both Horton and Roberts attempted to measure Ar. The mean of 
their values was 

A r = -0.0075 (12) 

In a later publication, 4 however, Roberts explains that the best agreement with experi- 
mental data was obtained by using Horton’s original theoretical value for the dissipation 
coefficient, 

C D = 0.022 (13) 


such that 


Ar = -.0059 


(14) 


Clearly, there are a number of contradictions among these results. For the time being, the 
values in Eqs. (8) and (12) are used in the bubble model. 

Assuming a constant Co and Hzz, Eqs. (8) and (9), and a linear pressure recovery 
from transition to reattachment, the energy integral equation can be integrated to obtain 
the momentum thickness at reattachment: 


(£2)* = (£2)7 






(15) 


Eliminating (62)/? between Eqs. (7) and (15), 

In order to implement this result numerically, Ur is decreased in small increments from the 
value of Ut- At each step, £2 is calculated and it is checked whether the segment joining 
T to R intersects the inviscid velocity distribution as shown in Fig. 1. The turbulent 
calculations start at sr, with the size of the first step equal to XB. This is the distance 
from sr to the first airfoil coordinate after reattachment. The drag coefficient is then 
determined from the Squire- Young formula. 
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Preliminary Results and Discussion 


Given all the assumptions discussed above, the modified Horton model cannot be 
expected to perform well in its present form. In addition, most of the assumptions imple- 
mented are based on measurements made at high Reynolds numbers, and are not neces- 
sarily correct at low Reynolds numbers. A good example of the weakness of this model is 
given by an analysis of the Eppler 387 airfoil, designed to operate at Re — 300,000. The 
aerodynamic characteristics predicted using the original Eppler analysis and those from 
the program with the bubble model included are compared in Fig. 2. Some points from the 
recent measurements of McGhee, et al. 7 are also included. The original Eppler analysis 
slightly underpredicts the drag, possibly because it does not account for the drag due to 
the bubble. Inclusion of the modified Horton model, as detailed above, however, likewise 
does not reproduce the experimental data. In fact, at the lower lift coefficients the drag 
predicted with the bubble model is lower than that predicted without it. In the mid-C* 
range, the prediction using the bubble model exhibits a large region of turbulent separated 
flow. Near a = 10° from zero lift, the length of the turbulent separated region upstream 
of the trailing edge decreases and the lift increases again. 

The low drag at low CV s is believed to be a consequence of assuming no growth in 
momentum thickness from separation to transition. This assumption appears to be valid 
for Re > 1,000,000 but, as shown by the data of Brendel and Mueller, not at low Reynolds 
numbers. The cause for the high drag and separated flow in the mid-Cf range is not clear. 
While the bubble length does not change significantly from its length at the lower angles of 
attack, the actual upper surface pressure levels may be noticeably lower than the inviscid 
values, especially upstream of the bubble. In other words, using the value of Us from the 
inviscid velocity distribution leads to a higher gradient across the bubble and the rest of 
the boundary layer than is present in reality; hence the deterioration in performance. As 
can be seen from the transition curves, the decreasing size of the separated flow length at 
the high Ci s is caused by the shift in bubble location from mid-chord to the leading edge. 

In order to better observe the details of the boundary-layer solution which lead to 
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the behavior just discussed, a new series of plots was generated. For each angle of attack 
of interest, the airfoil analyzed, the velocity distribution, the boundary-layer data, and 
the viscous analysis summary are all shown together on one page. This makes it possible 
to view all the relevant relationships simultaneously. Fig. 3a contains the results of the 
analysis at a — 1.52°, the low end of the polar. The upper surface exhibits a fairly long 
mid-chord bubble, while the lower surface a short leading-edge bubble. The boundary- 
layer data is plotted in the form used by Eppler as log 10 {i?e5 2 ) vs. #32, as well as in the 
four plots of H 12, 82, #32 > and c/ vs. s/c. Moreover, the results from the original Eppler 
method and those from the method with the bubble model are plotted together to aid in 
the comparison. To help in the interpretation of the data, two modifications were made to 
the original Eppler and Somers program. First, the actual separation and reattachment 
points are plotted in both the log and the linear plots and second, in the log plot, all 
the points for which H32 > 1-62 are plotted with H32 = 1-62 rather than being skipped 
altogether. 

At a = 2.6°, Fig. 3b, the character of the bubble on the upper surface has not 
changed, while the flow on the lower is all laminar. At a = 5°, Fig. 3c, it is clear that 
the model is not incrementing the momentum thickness correctly. In fact, the value at 
reattachment equals the value obtained by assuming transition at the laminar separation 
point. Incorporating an empirical correlation for the growth in 62 in the laminar part of 
the bubble should improve the results. At a — 8°, Fig. 3d, the bubble causes the flow 
to separate near the trailing edge. The higher drag coefficient is due to the separation 
as well as to an increase in 62 across the bubble which is higher than that through the 
same length of turbulent boundary layer. The dependence of the increment in 62 on 
angle of attack is related to the increasing steepness the inviscid pressure recovery as ct 
increases, the length of the bubble remaining approximately constant. At a = 10.48°, Fig. 
3e, this prediction trend is most noticeable. The fact that now the drag is even higher 
than the experimental values might be explained by the inability of this model to account 
for viscous/inviscid interaction. In fact, as long as the separated shear layer reattaches, 
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any adverse pressure gradient will drive the growth in boundary-layer properties across 
the bubble with no chance for them to “backfire” and mitigate the too severe pressure 
gradient. As shown in Fig. 23c in Ref. 7, however, the inviscid pressure distribution, even 
at a = 7.36°, has been modified. In particular, the rooftop, which is immediately upstream 
of the bubble, is lower, while the pressure at reattachment is closer to the inviscid value. 
On the other hand, from the observed data, a long bubble does not appear to be present, 
in the sense that the inviscid velocity distribution is not modified drastically. Therefore, 
in order to correctly capture the bubble drag, it may be sufficient to decrease the value of 
the velocity at separation obtained from the inviscid pressure distribution by means of an 
empirical factor. Such factor could be correlated to Reynolds number and angle of attack. 
This approach assumes that the history of the laminar boundary layer from stagnation to 
separation is not of primary importance. If such an assumption were to prove faulty, it 
might then be left as the only alternative to perform a displacement thickness iteration 
assuming transition at laminar separation, and then to apply a short-bubble model to the 
resulting “viscous” pressure distribution. At a = 10.49°, Fig. 3f, the laminar boundary 
layer separates upstream, near the leading edge, and a short bubble again forms. 

A further aspect of this model that will need to be treated correctly is the match 
to Eppler’s turbulent boundary-layer method at reattachment. Horton’s model, and the 
reattachment parameter A r, depend on the experimental values for (# 32 )/? = 1.51 and 
= 3.5. Eppler, however, uses a correlation for H 12 = # 12 (# 32 ) that yields H 12 = 
2.34 if H 32 — 1.51. This can be seen in all the plots described above. In other words, no 
allowance is made for that region, downstream of reattachment, in which the turbulent 
boundary layer is redeveloping and approaching a state for which Eppler’s correlations are 
valid. The importance of an accurate modelling of this region is not yet known. 
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Nomenclature 


c = Airfoil chord 

Cf = Skin-friction coefficient 

t\ = Laminar length of the bubble 

ti = Turbulent length of the bubble 

s = Streamwise coordinate from the stagnation point 

Cd = Sectional drag coefficient 

Cd = Dissipation coefficient 

Ci = Sectional lift coefficient 

C m = Sectional moment coefficient about c /4 

H 12 = Boundary layer shape factor, (Si/ 62) 

H32 = Boundary layer shape factor, (63/62) 

Ree, t = Transition Reynolds number, (Usti/u) 

Res 7 = Momentum thickness Reynolds number, (U62/v)s 

U = Velocity at the edge of the boundary layer 

Uoo = Free sream velocity 

XB = Distance from sr to first airfoil coordinate after sr 

XU = Laminar fraction within one boundary layer step 

a = Angle of attack relative to the zero-lift line 

£1 = Boundary layer displacement thickness 

62 = Boundary layer momentum thickness 

63 = Boundary layer kinetic energy thickness 

v = Kinematic viscosity of air 

A r = Horton’s reattachment parameter 

Subscripts: 

S = Laminar separation point 

R = Turbulent reattachment point 

T = Transition point 
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Fig. 1. Sketch of laminar separation bubble with important variables 
defined. 
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Aerodynamic characteristics for the Eppler 387 airfoil obtained 
with the original Eppler and Somers method compared with those 
obtained using the program incorporating the Horton short-bubble 
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Fig. 3d. a = 8.00 deg. from zero lift. 
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Fig. 3f. a * 10.49 deg. from zero lift. 








